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Abstract 

We present a new algorithm for reconstructing an exact algebraic num- 
ber from its approximate value using an improved parameterized integer 
relation construction method. Our result is almost consistent with the 
existence of error controlling on obtaining an exact rational number from 
its approximation. The algorithm is applicable for finding exact minimal 
polynomial by its approximate root. This also enables us to provide an 
efficient method of converting the rational approximation representation 
to the minimal polynomial representation, and devise a simple algorithm 
to factor multivariate polynomials with rational coefficients. 
Compared with other methods, this method has the numerical compu- 
tation advantage of high efficiency. The experimental results show that 
the method is more efficient than identify in Maple 11 for obtaining an 
exact algebraic number from its approximation. In this paper, we com- 
pletely implement how to obtain exact results by numerical approximate 
computations. 
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1 Introduction 



Symbolic computations are principally exact and stable. However, they have 
the disadvantage of intermediate expression swell. Numerical approximate com- 
putations can solve large and complex problems fast, whereas only give ap- 
proximate results. The growing demand for speed, accuracy and reliability in 
mathematical computing has accelerated the process of blurring the distinction 
between two areas of research that were previously quite separate. Therefore, 
algorithms that combine ideas from symbolic and numeric computations have 
been of increasing interest in recent two decades. Symbolic computations is 
for sake of speed by intermediate use of floating-point arithmetic. The work 
reported in [Jl [2 [3 |4l 13 [6] studies recovery of approximate value from numer- 
ical intermediate results. A somewhat related topic are algorithms that obtain 
the exact factorization of an exact input polynomial by use of floating point 
arithmetic in a practically efficient technique 7,8. In the meantime, symbolic 
methods are applied in the field of numerical computations for ill-conditioned 
problems [9l [10l [11] . The main goal of hybrid symbolic- numeric computation is 
to extend the domain of efficiently solvable problems. However, there is a gap 
between approximate computations and exact results [T2"]. 

We consider the following question: Suppose we are given an approximate 
root of an unknown polynomial with integral coefficients and a bound on the 
degree and size of the coefficients of the polynomial. Is it possible to infer 
the polynomial and its exact root? The question was raised by Manuel Blum 
in Theoretical Cryptography, and Jingzhong Zhang in Automated Reasoning, 
respectively. Kannan et al answered the question in [13] . However, their tech- 
nique is based on the Lenstra-Lenstra-Lovasz(LLL) lattice reduction algorithm, 
which is quite unstable in numerical computations. The function MinimalPoly- 
nomial in maple, which finds minimal polynomial for an approximate root, was 
implemented using the same technique. In this paper, we present a new algo- 
rithm for finding exact minimal polynomial and reconstructing the exact root 
by approximate value. Our algorithm is based on the improved parameterized 
integer relation construction algorithm, whose stability admits an efficient im- 
plementation with lower run times on average than the former algorithm, and 
can be used to prove that relation bounds obtained from computer runs using 
it are numerically accurate. The other function identify in maple , which finds 
a closed form for a decimal approximation of a number, was implemented using 
the integer relation construction algorithm. However, the choice of Digits of 
approximate value is fairly arbitrary [14] . In contrast, we fully analyze numer- 
ical behavior of an approximate to exact value and give how many Digits of 
approximate value, which can be obtained exact results. The work is regard 
as a further research in [15j . We solve the problem, which can be described as 
follows: 

Given approximate value a at arbitrary accuracy of an unknown algebraic 
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number, and we also know the degree of the algebraic number n and an upper 
bound N of its height on minimal polynomial in advance. The problem will be 
solved in two steps: First, we discuss how much control error e is, so that we 
can reconstruct the algebraic number a from its approximation a when it holds 
that \a — a\ < e. Of course, £ is a function in n and N. Second, we give an 
algorithm to compute the minimal polynomial of the algebraic number. 

We are able to extend our results with the same methods to devise a simple 
polynomial-time algorithm to factor multivariate polynomials with rational co- 
efficients, and provide a natural, efficient technique to the minimal polynomial 
representation. 

The rest of this paper is organized as follows. Section 2 illustrates the im- 
proved parameterized integer relation construction algorithm. In Section 3, we 
discuss how to recover a quadratic algebraic number and reconstruct minimal 
polynomial by approximation. Section 4 gives some experimental results. The 
final section concludes this paper. 



2 Preliminaries 

In this section, we first give some notations, and a brief introduction on integer 
relation problems. Then an improved parameterized integer relation construc- 
tion algorithm is also reviewed. 

2.1 Notations 

Throughout this paper, Z denotes the set of the integers, Q the set of the 
rationals, R the set of the reals, 0(R n ) the corresponding system of ordinary 
integers, U(n — 1,R) the group of unitary matrices over R, GL(n,0(M.)) the 
group of unimodular matrices with entries in the integers, co^B the i-th column 
of the matrix B. The ring of polynomials with integral coefficients will be denoted 
Z[A]. The content of a polynomial p(X) in Z[A] is the greatest common divisor 
of its coefficients. A polynomial in Z[X] is primitive if its content is 1. A 
polynomial p(X) has degree d if p(X) = ^2 i=0 PiX l with pd ^ 0. We write 
deg(p) = d. The length \p\ of p(X) — ^2h-qPiX 1 is the Euclidean length 
of the vector (po,Pi, • • • ,Pd)\ the height \p\oo of p(X) is the Loc-norm of the 
vector (p , .Pi, ■ ■ ■ ,Pd), so |p|oo = max <i<c2 \pi\. An algebraic number is a root of 
a polynomial with integral coefficients. The minimal polynomial of an algebraic 
number a is the irreducible polynomial in Z[X] satisfied by a. The minimal 
polynomial is unique up to units in Z. The degree and height of an algebraic 
number are the degree and height, respectively, of its minimal polynomial. 
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2.2 Integer relation algorithm 



There exists an integer relation amongst the numbers x\,x%, ■ ■ ■ ,x n if there are 
integers 01,02,- •• ,a n , not all zero, such that X)"=i °>i x i = 0- For the vector 
x = [xi, X2, • • • , x n ] T , the nonzero vector a € Z™ is an integer relation for x if 
a • x = 0. Here are some useful definitions and theorems [TBI 1 1 Tj : 

Definition 1 (M x ) Assume x — [xi,X2, ■ ■ • ,x n ] T £ E™ has norm \x\=l. De- 
fine x 1 - to be the set of all vectors in M. n orthogonal to x. Let 0(R n ) fl x^be 
the discrete lattice of integral relations for x. Define M x > to be the smallest 
norm of any relation for x in this lattice. 

Definition 2 (H x ) Assume x = [xi,X2,- ■ ■ ,x n ] T S M. n has norm \x\=l. Fur- 
thermore, suppose that no coordinate entry of x is zero, i.e., Xj ^ for 1 < j < 
n(otherwise x has an immediate and obvious integral relation). For 1 < j < n 
define the partial sums 



Given such a unit vector x, define the n x (n — 1) lower trapezoidal matrix 
H x = (hi,j) by 



[-XiXj/(sjS j+1 ) if 1 < j < i < n. 
Note that hij is scale invariant. 

Definition 3 (Modified Hermite reduction) Let H be a lower trapezoidal matrix, 
with hij = if j > i and hjj 7^ 0. Set D = I n , define the matrix D = 
[dij) G GL(n, 0(i?)) recursively as follows: For i from 2 to n, and for j from 
i-1 to l(step-l), set q = nint(hi t j/hjj); then for k from 1 to j replace hi^k by 
hi t k — qhj t k, and for k from 1 to n replace d^u — qdj,k> where the function nint 
denotes a nearest integer function, e.g., nint(t)=\t + 1/2J. 

Theorem 1 Let x =/= € M". Suppose that for any relation m of x and for 
any matrix A £ GL(n,0(M)) there exists a unitary matrix Q£ U(n-l) such that 
H = AH X Q is lower trapezoidal and all of the diagonal elements of H satisfy 
hj j ^ 0. Then 




j<k<n 




ifl<i<j<n — 1, 
tfl<i = j<n-l, 



1 . 1 



max l<j<n-l|/y,i| 1<J<«-1 \hjj\ 
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Remark 1 The inequality of Theorem [7] offers an increasing lower bound on 
the size of any possible relation. Theorem\j] can be used with any algorithm that 
produces GL(n,0(M)) matrices. Any GL(n,Q(M)) matrix A whatsoever can be 
put into Theorem^ 

Theorem 2 Assume real numbers, n > 2, r > 1, 7 > \J 4/3, and that 7^ 
x G R" /ias O(K) integer relations. Let M x be the least norm of relations for x. 
Then PSLQ{t) will find some integer relation for x in no more than 

\2j logr 

iterations. 

Theorem 3 Let M x be the smallest possible norm of any relation for x. Let m 
be any relation found by PSLQ(t). For all 7 > ^4/3 for real vectors 

\m\ < j n - 2 M x . 

Remark 2 For n=2, Theorem^ proves that any relation O^rnG 0(M 2 ) found 
has norm \m\ — M x . In other words, PSLQ(t) finds a shortest relation. For 
real numbers this corresponds to the case of the Euclidean algorithm. 

Based on the theorems as above, and if there exists a known error controlling 
e, then an algorithm for obtaining the integer relation was designed as follows: 

Algorithm 1 Parameterized integer relation construction algorithm 
Input: a vector x, and an error control e > 0; 
Output: an integer relation m. 

Step 1: Set i := 1, r := 2/\/3, and unitize the vector x to x; 
Step 2: Set H x by definition^ 

Step 3: Produce matrix D £ GL(n, O(R)) using modified Hermite Re- 
duction; 

Step 4: Set x := x ■ D~ l ,H := D ■ H, A := D ■ A, B := B ■ D -1 , 
case 1: if xj — 0, then m := coljB, goto Step 11; 
case 2: if 'ha < e, then m := col n _\B, goto Step 11; 

Step 5: i := i + 1; 

Step 6: Choose an integer r, such that T r |/i rjr | > T 3 |/ij j|, for all 1 < 
3 <n - 1/ 
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Step 7: Define a :— h r>r , (3 := h r+ i_ r , A := /i r +i,r+i, o~ '■= \J fl 2 + A 2 ; 

Step 8: Change h r to /i r +i, and define the permutation matrix R; 

Step 9: Set x := x ■ R, H := R- H, A := R- A, B := B ■ R, if i=n-l, 
then goto Step 4; 

Step 10: Define Q := (q hj ) G U(n - 1, R), H := H ■ Q, goto Step 4; 
Step 11: return m. 

By algorithm!]] w e can find the integer relation m of the vector x = (f , a, a 2 , ■ ■ ■ , ct n ). 
So, we get a nonzero polynomial of degree n, i.e., 

G{x) =m- {l,x,x 2 ,--- ,x n ) T . (1) 

Our main task is to show that polynomial |T]) is uniquely determined under 
assumptions, and discuss the controlling error e in algorithm [T] in the next 
section. 



3 Reconstructing minimal polynomial from its 
approximation 



In this section, we will solve such a problem: For a given floating number a, 
which is an approximation of unknown algebraic number, how do we obtain 
the exact value? Without loss of generality, we first consider the recovering 
quadratic algebraic number from its approximate value, and then generalize the 
results to the case of algebraic number of high degree. At first, we have some 
lemmas as follows: 

Lemma 1 Let f = ^ZILo a i xl ^ %[x] be a polynomial of degree n > 0, and let 
e = maxi<i<„ \a l — &i\ for the rest of this paper, where di for 1 < i < n are the 
rational approximations to the powers a % of algebraic number a, and do = 1. 
Then 

\f(a)-f(d)\ <e-n -l/loo- (2) 

Proof: Since f(a)-f(d) = £" =0 a^af-di), we get |/(a)-/(5)| = |E"=i a i( ai - 
di)\, and then 

n n n 

| a t (a l - di)\ < N ' - on)\ < \ai\ ■ e < n ■ l/U • e. 

i—l i—1 i—1 

So, the lemma is finished. 
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Lemma 2 Let h and g be nonzero polynomials in Z[x] of degree n and m, 
respectively, and let a E K be a zero of h with \a\ < 1. If h is irreducible and 
g(a) ^ 0, then 

|g(a)|>n- 1 -H~ m -M 1 - n . (3) 
Proof: See Proposition(1.6) of[13]. If \a\ > 1, a simple transform of it does. 

Corollary 1 Let h and g be nonzero polynomials in Z[x] of degrees n and m, 
respectively, and let a E R be a zero of h with \a\ < 1. //ft. is irreducible and 
g{a) ^ 0, i/ien 

|.g(a)| > n- 1 • („ + 1)^ • (m + l) 1 ^ • \h\^ ■ \g\^ n . (4) 

Proof: First notice that |/| 2 < (n + \) ■ \ f\ 2 OQ holds for any polynomial / of 
degree at most n > 0, so |/| < \/n + 1 • |/|oo- From LemmaUwe get 

| 5 (a)| > n- 1 • (n + 1)"* • (m + l) 1 ^ • \h\~™ • |. 9 |L~"- 
So, the corollary is finished. 

Theorem 4 Lei cm approximate value a belong to an unknown algebraic num- 
ber a of degree n > 0. Assume that the existence of the polynomial G{x) = 
Y^i=o a i xl > where a n ^ 0. Suppose n and upper bound N on the degree and 
height of minimal polynomial g(x) on the algebraic number a are known, re- 
spectively. If 

\G(a)\ < n- 1 ■ (n + 1)-"+* • \G\ x n ■ N^ n - n ■ e ■ |G|«,, 

then 

\G(a)\ <n- 1 -(n+l)- n+ '-|G|-"-iV 1 - n . 

Proof: From lemma [TJ we notice that \G(a) — G(a)\ < e ■ n ■ |G|oo, and 

\G(a) - G{a)\ > \G(a)\ - \G(a)\, so \G(a)\ < \G(a)\ + n ■ e ■ (G^. 
From the assumption of the theorem, since 

|G(S)| < n" 1 ■ (n + l)- n+ i • \G\ x n ■ N^ n - n ■ e ■ |GU (5) 
So, the theory is finished. 

Corollary 2 If \G(a)\ < n" 1 • (n + l)- n+ i • |G|^ n • N 1 ^, where G(x) is con- 
structed by the parameterized integer relation construction algorithm as above, 
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the upper bound N on the height of its minimal polynomial g(x) on an algebraic 
number a are known. Then 

G(a) = 0, (6) 

and the primitive part of polynomial G(x) is the minimal polynomial of algebraic 
number a. 

Proof: Proof is given by contradiction. According to Lemma [H suppose on the 
contrary that G(a) ^ 0, then 

\G(a)\>n- l -(n+l)- n+ ? ■ \G\^ ■ N^ n . 

From theory |U we get 

\G(a)\<n- x .(n + l)- n+ ? • |G| M " • TV 1 "". 

So, G(a) = 0. Since algebraic number a is degree n > 0, then the primitive 
polynomial of G{x) denotes pp(G(x)), hence pp(G(x)) is just irreducible and 
equal to g{x). Of course, it is unique. 

So, the corollary is finished. 

3.1 Recovering quadratic algebraic number from approx- 
imate value 

For simplicity, we discuss how to obtain quadratic algebraic number from its 
approximation using integer relation algorithm. Let a be the approximate value, 
considering the vector v = (l,a,a 2 ). Our goal is to find a vector w which has 
all integer entries such that the dot product of v and w is less than a lower 
bound, which is obtained and we are able to get the size of the neighborhood is 
l/(12v / 3-/V 4 ) from theorem[U The following theorem answers the basic questions 
of this approach. 

Theorem 5 Let a be an approximate value belonging to an unknown quadratic 
algebraic number a, if 

e = \a-a\ < l/(12 v / 3A r4 ), (7) 

where N is the upper bound on the height of its minimal polynomial. Then 
G(a) — 0, and the primitive part of G(x) is its minimal polynomial, where 
G{x) = X)i=o aixl * s constructed using integer relation algorithm as above. 

Proof: From theorem [I] and corollary O if and only if 

\G[a)\ <n- 1 .(n+l)-"+^.|G| TC "-7V 1 -», 
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Therefore, G{a) — 0. Under the assumption of the theorem we get n — 2 and 
|G(a)| > 0, hence it is obvious that inequality ([5]) holds, i.e., 

< n- 1 ■ (n + • |G|-" • AT 1 "" -n-e- \G\oo. (8) 

So, Solving inequality ([8]) yields 

e<l/(l2V3\G\lN). 

From theorem [3j and a is a quadratic algebraic number, so |G|oo is just equal 
to N. 

So, the theory is finished. 

This theorem leads to the following algorithm for recovering the quadratic 
algebraic number of a: 

Algorithm 2 Recovering quadratic algebraic number algorithm 

Input: an floating number (a, N) belonging to an unknown quadratic algebraic 

number a, i.e., satisfying |?]]. 

Output: an quadratic algebraic number a. 

Step 1: Construct the vector v; 

Step 2: Compute e satisfying 

Step 3: Call algorithm^ to find an integer relation w for v; 

Step 4 : Obtain w(x) the corresponding polynomial; 

Step 5: Let g{x) be the primitive part of w{x); 

Step 6: Solve the equation g(x) = and choose the corresponding 
algebraic number to a; 

Step 7: return a. 

Theorem 6 Algorithm^ works correctly as specified and uses OflogN) binary 
bit operations, where N is the upper bound of height on its minimal polynomial. 

Proof: Correctness follows from theorem [5j The cost of the algorithm is 
0(logN) binary bit operations obviously. 
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3.2 Obtaining minimal polynomial of high degree 

If a is a real number, then by dcfinitioin a is algebraic exactly if, for some n, 
the vector 

(l,a,« 2 ,"-,a n ) (9) 

has an integer relation. The integer coefficient polynomial of lowest degree, 
whose root a is, is determined uniquely up to a constant multiple; it is called 
the minimalpolynomial for a. Integer relation algorithm can be employed to 
search for minimal polynomial in a straightforward way by simply feeding them 
the vector ([9]) as their input. Let a be an approximate value belonging to an 
unknown algebraic number a, considering the vector v — (1, a, a 2 , ■ ■ ■ , a"), how 
to obtain the exact minimal polynomial from its approximate value? We have 
the same technique answer to the question from the following theorem. 

Theorem 7 Let a be an approximate value belonging to an unknown algebraic 
number a of degree n > 0, if 

e=\a-&\<l/(n 2 {n + l) n -iN 2n ), (10) 

where N is the upper bound on the height of its minimal polynomial. Then 
G(a) — 0, and the primitive part of G(x) is its minimal polynomial, where 
G(x) = X^r=o aixl * s constructed using the parameterized integer relation con- 
struction algorithm as above. 

Proof: The proof can be given similarly to that in theorem [5j 

It is easiest to appreciate the theorem by seeing how it justifies the following 
algorithm for obtaining minimal polynomials from its approximation: 

Algorithm 3 Obtaining minimal polynomial algorithm 

Input: an floating number (a, n, N) belong to an unknown algebraic number a, 
i.e., satisfying mU\) . 

Output: g(x), the minimal polynomial of a. 

Step 1: Construct the vector v; 

Step 2: Compute e satisfying \l(fy : 

Step 3: Call algorithm^ to find an integer relation w for v; 

Step 4 : Obtain w(x) the corresponding polynomial; 

Step 5: Let g{x) be the primitive part of w(x); 

Step 6: return g(x). 
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Theorem 8 Algorithm\3\works correctly as specified and uses 0(n(logn+logN) ) 
binary bit operations, where n and N are the degree and height of its minimal 
polynomial, respectively. 

Proof: Correctness follows from theorem [Jj The cost of the algorithm is 
0(n(logn + logN)) binary bit operations obviously. 

The method of obtaining minimal polynomials from an approximate value can 
be extended to the set of complex numbers and many applications in computer 
algebra and science. 

This yields a simple factorization algorithm for multivariate polynomials with 
rational coefficients: We can reduce a multivariate polynomial to a bivariate 
polynomial using the Hilbert irreducibility theorem, the basic idea was described 
in [5], and then convert a bivariate polynomial to a univariate polynomial by 
substituting a transcendental number in [18] or an algebraic number of high 
degree for a variate in [H]. It can find the bivariate polynomial's factors, from 
which the factors of the original multivariate polynomial can be recovered using 
Hensel lifting. After this substitution we can get an approximate root of the 
univariate polynomial and use our algorithm to find the irreducible polynomial 
satisfied by the exact root, which must then be a factor of the given polynomial. 
This is repeated until the factors are found. 

The other yields an efficient method of converting the rational approximation 
representation to the minimal polynomial representation. The traditional repre- 
sentation of algebraic numbers is by their minimal polynomials jTS} [5j]] [5TJ [55] . 
We now propose an efficient method to the minimal polynomial representation, 
which only needs an approximate value, degree and height of its minimal poly- 
nomial, i.e., an ordered triple < a,n,N > instead of an algebraic number a, 
where a is its approximate value, and n and N are the degree and height of its 
minimal polynomial, respectively, denotes < a >—< a, n, N >. It is not hard to 
see the computations in the representation can be changed to computations in 
the other without loss of efficiency, the rational approximation method is closer 
to the intuitive notion of computation. 



4 Experimental Results 

Our algorithms are implemented in Maple. The following examples run in the 
platform of Maple 11 and PIV3.0G,512M RAW. The first three examples illumi- 
nate how to obtain exact quadratic algebraic number and minimal polynomials. 
Example 4 tests our algorithm for factoring primitive polynomials with integral 
coefficients. 

Example 1 Let a be an unknown quadratic algebraic number. We only 
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know an upper bound of height on its minimal polynomial N — 47. According 
to theorem [5j compute quadratic algebraic number a as follows: First obtain 
control error e = 1/(12 * y/3 * TV 4 ) = 1/(1807729447692 * V3) « 1.0 x 10~ 8 . 
And then assume that we use some numerical method to get an approximation 
a = 11.937253933, such that \a — a\ < e. Calling algorithm [2] yields as follows: 
Its minimal polynomial is g(x) — x 2 — 8 * x — 47. So, we can obtain the 
corresponding quadratic algebraic number a = 4 + 3\/7. 

Remark 3 The function identify in maple 11 needs Digits^ 13, whereas our 
algorithm only needs 9 digits. 

Example 2 For obtaining exact minimal polynomials from approximate root 
a, we only know degree n = 3 and height N = 17 of its minimal polynomial. 
According to theorem just as do in Example 1: First get the error e — 
l/{n 2 (n+l) n -^N 2n ) = 1/6951619872 « 1.4xl0~ 10 . Assume that we use some 
numerical method to get an approximation a = 16.808034642702, such that 
\a — a | < e. Calling algorithm [3] yields as follows: Its minimal polynomial is 
g{x) = x 3 - 17 * x 2 + 4 * x - 13. 

Example 3 Let a known floating number a belonging to some algebraic 
number a of degree n = 4, where a — 3.14626436994198, we also know an upper 
bound of height N = 10 on its minimal polynomial. According to theorem [71 
we can get the error e = l/{n 2 {n + l) n ~^ N 2n ) = l/(4 2 * 5^ * 10 8 ) «2.2xl0" 12 . 
Calling algorithm^ if only the floating number a, such that \a — a\ < e, then 
we can get its minimal polynomial g(x) = x 4 — 10 * x 2 + 1. So, the exact 
algebraic number a is able to denote < a >=< 3.14626436994198,4, 10 >, i.e., 
< V2 + VS >=< 3.14626436994198, 4, 10 >. 

Example 4 This example is an application in factoring primitive polynomials 
over integral coefficients. For the conveniency of display in the paper, we choose 
a very simple polynomial as follows: 

p = 3x 9 - 9x s + 3x 7 + 6x 5 - 27x 4 + 21x 3 + 30x 2 - 2lx + 3 

We want to factor the polynomial p via reconstruction of minimal polynomials 
over the integers. First, we transform p to a primitive polynomial as follows: 

p = x 9 - 3x 8 + x 7 + 2x 5 - 9x 4 + 7x 3 + 10x 2 -7x+l, 

We see the upper bound of coefficients on polynomial p is 10, which has relation 
with an upper bound of coefficients of the factors on the primitive polynomial 
p by Landau-Mignotte bound [23]. Taking N = 5, n = 2 yields e = l/(2 2 * (2 + 
l) 2 -3 * 5 4 ) = 1/(7500 * V$) w 8.0 x 10~ 5 . Then we compute the approximate 
root on x. With Maple we get via [fsolve(p = 0,a;)]: 

S = [2.618033989, 1.250523220, -.9223475138, .3819660113, .2192284350] 
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According to theorem[71 let a — 2.618033989 be an approximate value belonging 
to some quadratic algebraic number a, calling algorithm [3] yields as follows: 

Pi = x 2 — 3 * x + 1. 

And then we use the polynomial division to get 

P2 = x 7 + 2 * x 3 — 3 * x 2 — 4 * x + 1. 

Based on the Eisenstein's Criterion [24] , the P2 is irreducible in Z[A]. So, the 
Pi and pi are the factors of primitive polynomial p. 



5 Conclusion 

In this paper, we have presented a new method for obtaining exact results by 
numerical approximate computations. The key technique of our method is based 
on an improved parameterized integer relation construction algorithm, which is 
able to find an exact relation by the accuracy control e in formula (|10p is an 
exponential function in degree and height of its minimal polynomial. The re- 
sult is almost consistent with the existence of error controlling on obtaining an 
exact rational number from its approximation in [15] . Using our algorithm, we 
have succeed in factoring multivariate polynomials with rational coefficients and 
providing an efficient method of converting the rational approximation represen- 
tation to the minimal polynomial representation. Our method can be applied in 
many aspects, such as proving inequality statements and equality statements, 
and computing resultants, etc.. Thus we can take fully advantage of approxi- 
mate methods to solve larger scale symbolic computation problems. 
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